Bond order from disorder in the planar pyrochlore magnet 
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We study magnetic order in the Heisenberg antiferromagnet on the checkerboard lattice, a two- 
dimensional version of the pyrochlore network with strong geometric frustration. By employing the 
semiclassical (1/S) expansion we find that quantum fluctuations of spins induce a long-range order 
that breaks the four-fold rotational symmetry of the lattice. The ordered phase is a valence-bond 
crystal. We discuss similarities and differences with the extreme quantum case S — 1/2 and find a 
useful phenomenology to describe the bond-ordered phases. 



I. INTRODUCTION 

Frustrated magnets have attracted the attention of 
theorists for several decades, beginning with the study 
of the Ising antiferromagnet on the triangular latticed 
More recently, new families of frustrated magnetic com- 
pounds have become available for experimental studies 
reviving the interest in their properties'^ By its very 
nature, a frustrated system has an extremely large clas- 
sical degeneracy of the ground state. This degeneracy 
is very effective in suppressing classical spin order ^ thus 
providing a route to non-Neel (quantum) ground states 
even for higher-dimensional systems. The nature of such 
ground states is far from obvious: the aforementioned 
degeneracy allows for a variety of unusual vacua. Among 
the possibilities are bond-ordered states, in which pair av- 
erages (Si-Sj), rather than spins (Si) themselves, form a 
periodic patternjrSil and spin liquids that break no lattice 
symmetry but are distinguished by the unusual quantum 
numbers and statistics of their excitations^ 

One of the most intensively studied frustrated systems 
is the Heisenberg antiferromagnet on the pyrochlore lat- 
tice. It has many experimental realizations that show 
rather remarkable magnetic behavior. For example, the 
spinel ZnCr204 is the first frustrated magnet in which 
zero modes — spin waves connecting degenerate ground 
states — have been observed by neutron scattering^ Not 
long ago, it has been showniSiii that a coupling between 
spins and lattice vibrations leads to a spin-Peierls phase 
transition in this manifestly three-dimensional spin sys- 
tem. This effect, also observed^ in ZnCr204, is classical 
in nature in the sense that it is not parametrically small 
in 1/S. Therefore it is expected to dominate the more 
subtle quantum effects for large values of spin. 

Effects of frustration in quantum pyrochlore antiferro- 
magnets, particularly in the limit of a small spin 5 1 , are 



drawing quite a bit of interest. Finding answers in this 
case may provide clues to the unusual behavior of un- 
derdoped cuprate superconductors, where frustration of 
the spin system is achieved through the motion of doped 
charges. Although replacing dynamic frustration with ge- 
ometric one somewhat simplifies the problem, it is still far 
from trivial. An exact solution for S = 1/2 is not avail- 
able and is not expected in the immediate future. Numer- 
ical diagonalizations are hampered by the quick growth 
of the Hilbert space with the lattice size in three dimen- 
sions. Several research groups are attacking the prob- 
lem from various solvable limits: large 5^ large 
and weakly-coupled spin clusters. 15 i 16 i 17 i 18 i 19 Because it 
is not even obvious that extrapolations from these limits 
will lead to a converging answer, it seems highly desirable 
to test these approaches on a similar problem for which 
numerical answers are available. 

Most recently, a two-dimensional version of the py- 
rochlore network, the checkerboard lattice^ (also known 
as the planar pyrochlore and the square lattice with 
crossings), has become a focus of analytical2i*22i2^ and 
numerical22i2£iSii2£ studies. Lower dimensionality of 
this system makes it an easier target for numerical ap- 
proaches; at the same time, it has the local coordination 
of the pyrochlore lattice: magnetic bonds form a network 
of corner-sharing tetrahedra with spins at the vertices. It 
is therefore hoped that studies of the Heisenberg antifer- 
romagnet on the checkerboard lattice can shed light on 
the behavior of its three-dimensional analog. 

The planar pyrochlore lattice does differ from the py- 
rochlore proper in one fundamental aspect: not all bonds 
of its tetrahedra are equivalent — because no symmetry of 
the lattice turns first neighbors (horizontal and vertical 
bonds) into second (diagonal) . Even if the corresponding 
exchange couplings J\ and J 2 are set equal, spin corre- 
lations between first and second neighbors tend to be 
different, as evidenced by both analytical and numerical 
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FIG. 1: Heisenberg antiferromagnet on the checkerboard lat- 
tice. Shown is a generic classical ground state for a model 
in which the second-neighbor coupling J2 that exceeds the 
nearest-neighbor coupling J\ . At the level of the classical ap- 
proximation, the Neel order parameters of individual chains 
fii are uncorrelated. The direction is defined as that of 
spin Si at the left edge. 



results. A lack of such symmetry compels one to look at 
the general case with J x ^ J 2 . 

In this work, we study the checkerboard antiferromag- 
net in the limit of large spin 5, which allows for a sys- 
tematic perturbation theory in powers of 1/5. 28 We then 
compare our answers to the available numerical results 
for the opposite limit, 5 — 1/2, and find a simple phe- 
nomenology that describes both rather well. Our under- 
standing of the checkerboard antiferromagnet has been 
greatly helped by three recent ideas: (a) Henley's 'gauge 
symmetry' relating degenerate collinear ground states in 
a frustrated magnet ^ (b) casting of the problem in terms 
of bond — rather than spin — variables, which were intro- 
duced a decade ago by Harris et al.f (c) a realization 
of the potential significance of the q = 3 Potts model to 
bond-ordered states on the pyrochlore and checkerboard. 

Because this paper is long and technical, the reader 
may find it helpful to peruse an informal introduc- 
tion to the subject written by one of us~ That ar- 
ticle states our reasons to pursue the large- 5 route to 
strongly frustrated quantum magnets, explains the chal- 
lenges of that approach and points out ways to overcome 
them. It also contrasts the results obtained for differ- 
ent two-dimensional analogs of the pyrochlore lattice (the 
checkerboard is one of them) . 

The ground state of a classical (5 — 00) Heisenberg 
magnet is found by minimizing its energy 

E = Jij s * ' s 3 = ®(S 2 ) (!) 

(ij) 



with respect to classical spin variables S». Exchange cou- 
pling is Ji on horizontal and vertical bonds, and J 2 along 
diagonals. For weaker diagonal bonds, J2 < J\, classical 
energy minimization gives a unique ground state (modulo 
a global rotation of all spins). The ground state, shown 
in Fig. |2fa), is the same as that of the simple square 
lattice (J 2 = 0). 

In the region J 2 > Ji, the classical ground state be- 
comes continuously degenerate. For stronger diagonal 
bonds, J 2 > Ji, the system can be viewed as a collec- 
tion of criss-crossing chains running along diagonals of 
the lattice (Fig. QJ. Let us choose north-east to be our 
positive ^-direction, and north-west to be 77. Classically, 
each chain has perfect Neel order at zero temperature, 
however, directions of staggered magnetizations n.; of dif- 
ferent chains are completely independent at the classical 
level. In an L x L lattice with periodic boundary con- 
ditions classical ground states can be parametrized by L 
unit vectors n^. The classical degeneracy increases even 
further in the case of equal exchanges J 2 = JimQ£L 

The first-order (in 1/5) correction to the classical en- 
ergy comes from the zero-point quantum fluctuations of 
spin waves, 

£ 1 = 5>K|/2 = 0(5), (2) 

k 

where {ujk} are eigenfrequencies of classical spin waves 
about a given ground state obtained from the equations 
of motion 

//S, )J./,,S, • S,. (3) 

3 

It has been established previously that, quite generally, 
quantum fluctuations select ground states with collinear 
spins (assuming such classical ground states exist) 1 32 ' 33 

For a fixed global direction n, the problem thus reduces 
to a minimization of the zero-point energy J2J) over a dis- 
crete set of collinear Neel states. Thus selected ground 
states can be characterized in the thermodynamic limit 
with the aid of some order parameters. In addition to 
violating the spin- rotation symmetry 0(3), these ground 
states can also break some discrete lattice symmetries. 
For instance, when the ground states are not symmetric 
under 7r/2 rotations of the plane, one expects an order pa- 
rameter with the structure Z2 x S 2 . For a given direction 
of the Neel vector, there should then be two degenerate 
ground states. (See, e.g., the work by Chandra et alM- 
on the square lattice with a large second-neighbor cou- 
pling.) Contrary to these expectations, we find that for 
J 2 > J\ the ground state is fourfold degenerate with an 
order parameter Z 2 x Z 2 x S 2 . The extra degeneracy 
turns out to be related to a gauge-like symmetry (due to 
HenleyJ^,) that exists at the order 1/5 in the semiclassi- 
cal expansion. This dynamical symmetry is responsible 
for an even higher degeneracy of the ground state at the 
most frustrated point J\ = J 2 . In that case, the Neel 
order is destroyed and the order parameter is reduced 
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to Zi x Z 2 x S 2 / Z2 (two Ising orders and a director). 
At Ji = J2 the planar pyrochlore is a valence-bond solid 
with two independent bond orders and a nematic spin 
order. 

The paper is organized in the following way. In most 
of it ( Sections ITTI through E} we study the lowest-order — 
0(1/S) — quantum corrections to the degenerate classi- 
cal limit. We first explore the case where the second- 
neighbor coupling J 2 dominates and the system can 
be viewed as a set of weakly coupled antiferromagnetic 
chains. The interchain coupling J\ is frustrated and 
has no effect at the classical level. We show in Sec. [n] 
that, in line with the standard arguments^S quan- 
tum fluctuations favor collinear spin states. By using 
J1/J2 as a small parameter, we derive an effective in- 
teraction between the chains generated by quantum fluc- 
tuations. This potential is minimized by four distinct 
classical states. In Sec. lIIII we prove that this degeneracy 
remains intact for all J\ < Ji and trace its origins to Hen- 
ley's gauge symmetry. Sec. IIVI presents our findings in 
the strongly frustrated case of equal exchange couplings. 
This time, a much larger degeneracy of the ground state 
kills the Neel order (replacing it with a nematic order), 
but the bond order survives. A summary of the large-5 
results is given in Sec.[V] Finally, in Sec. I VII we explore 
the connection of these large- S results to the S — 1/2 
phase diagram obtained in numerical studies and specu- 
late on a phenomenology of the bond order in the S =1/2 
case. 





(a) 
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II. WEAKLY COUPLED CHAINS: Ji < J 2 

For Ji = 0, the magnet is reduced to a collection 
of independent antiferromagnetic chains running along 
the diagonals (Fig. 0). In a classical ground state, each 
chain has perfect antiferromagnetic order. For this rea- 
son, there is no coupling between intersecting chains at 
the classical level — even in the presence of a finite inter- 
chain coupling J 1 < J2. 

Quantum fluctuations disrupt the perfect Neel align- 
ment of adjacent spins and thus enable the chains to in- 
teract. For a weak interchain coupling J\ <C J2, one can 
use a systematic perturbation theory in Ji/J 2 (staying 
at the same order in 1/S) developed by Shenderi^ 



A. Effective interactions between chains 

To the lowest nontrivial order in J1/J2, the interchain 
coupling generates a potential that selects collinear spin 
configurations. For Neel magnetizations of individual 
chains n m , 



2tt J 2 



E 

crossings 



(n m -n„) 2 



(4) 



This interaction — coupling any two chains intersecting 
over a tetrahedron — is minimized when all staggered 



FIG. 2: (a) The classical ground state for Ji > J2. (b) 
through (e) are the four classical ground states at the order 
1/5 for Ji < J2. Frustrated bonds (those with two parallel 
spins) are shown in color dashed lines. Bottom figures show 
the lattice of tetrahedra; the primary colors (red, green, and 
blue) encode the location of frustrated bonds. See Fig. [J| for 
more details. 




FIG. 3: Bond variables f = {fi,fi} for a single tetrahedron 
in a ground state Si = 0. See Eq. 12UH for a definition 

of f. (a) Classical spins, (b) S = 1/2. In both cases, primary 
colors denote the pure states with two frustrated bonds (par- 
allel classical spins or a spin triplet). Secondary colors mark 
the pure states with two satisfied bonds (antiparallel spins 
or a spin singlet). The triangle and the circle delineate the 
domains of attainable values. 
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magnetizations point along a common axis, n m = ±n. 
The tendency of spins to align is the standard outcome 
of the order- from-disorder scenario i^SiS See Appendix for 
a derivation of Eq. 

Collinear Neel states can be characterized by a col- 
lection of Ising variables: s m = ±1 for chains running 
along the direction £ and t m — ±1 for chains along rj, 
see Fig. ^ Given a global orientation of spins n, these 
numbers determine the staggered magnetizations of in- 
dividual chains, n m = s m n or t m n, respectively. To the 
fourth order in J% (but still to the first order in 1/5) we 
obtain the correction to the classical energy of the ground 
state 



2JfS 
Jo 3 



2kS 



(5) 



2 i.— l ; — 



k=l 1=1 m n 



2l<>m<'m+2l 



B2k-l,2l-l S m S m +2H t n t n +2l-l) i 



where the positive coefficients Ai and Bki are computed 
in the Appendix. Note that there is no pairwise interac- 
tion between adjacent parallel chains (nor, for that mat- 
ter, between any parallel chains an odd distance apart). 
As a result of that, the effective interaction (J5J is invari- 
ant under the transformation 

Sm >-> {-iTSra, t n I * (-l)X, (6) 

which flips the spins on every other diagonal chain. 



on classical chains. On the other, they create interchain 
coupling and make the system two-dimensional thereby 
making spin order more likely. Which tendency wins? 

To answer this question one can compute the expec- 
tation values of local magnetization (S r ) in one of the 
ground states [Fig. |2fb) through (e)] . The classical value 
(0, 0, ±5) is reduced by quantum fluctuations of spins. 
A naive evaluation of this quantity at the first order 
in 1/5 gives a divergent negative correction suggesting 
that the order is destroyed. This, however, is an artefact 
of a low-order approximation. The magnon spectrum is 
given by the frequencies of classical spin waves, which 
know nothing about the interchain coupling. As a resut, 
the magnon spectrum contains lines of zero modes along 
the diagonal directions (pj = 0,7r and — 0,7r). The 
abundance of soft modes leads to a destruction of the 
long-range spin order. 

At the next level of approximation, C(l/5 2 ), magnon 
interactions modify the spin-wave spectrum lifting the 
zero modes to finite frequencies (with the exception of 
isolated points in the Brillouin zone).-- The spin exci- 
tations become two-dimensional and the infrared diver- 
gence of the correction to local magnetization is removed. 
Long-range order can survive. 

Alternatively, the hardening of the spin-wave spec- 
trum (specifically of the zero modes) can be evaluated 
already at the order 0(1/ S) by adding to the Heisenberg 
Hamiltonian a phcnomcnological biquadratic exchange 



ter 
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B. Ground states 



(10) 



The largest term in Eq. © is the interaction of 4 chains 
intersecting around an empty plaquette: 



2(7^5 v -v v -v 

2-^1 i' 1 SmSm + 1 tnt n +l 

J 2 ^ ' ' 



(7) 



(see Appendix for details). It is minimized by ground 
states of two distinct kinds: 

s m = Sm+i, t n = t n+ i, Figl^b) and (c), (8) 
s m = -«m+i, t n = -tn+i, Figl^d) and (e). (9) 

Curiously, the two types of ground states JHJ and © are 
related to each other by the staggering transformation 
©, rather than by any symmetry of the lattice. They 
remain degenerate even upon the inclusion of all two- 
and four-chain interactions in Eq. (JSJ . The origin of this 
dynamical symmetry will be discussed in Sec. IIIII 



C. Long-range spin order 

It is natural to ask whether there is a long-range spin 
order in the system. Quantum effects turn out to be 
rather subtle in this case. On the one hand, quantum 
fluctuations tend to destroy the long-range order found 



that couples spins at the intersections of chains. The 
strength of this interaction K = 0(Jf / J2S 3 ) is chosen 
so as to mimick, at the classical level, the energy of zero- 
point fluctuations. The latter is given, to the lowest order 
in Ji, by Eq. I|A.24|) . The spin- wave spectrum is then 
computed from the classical equations of motion for the 
spins. 

By performing a calculation along these lines (to 
be reported elsewhere 3 ^) we find that indeed the zero 
modes acquire energies of order Ji\f~S. The renormal- 
ized magnon frequency vanishes at the points 



(P£,Pr,) = (0,0) and 



(11) 



only, as required by the Goldstone theorem. The average 
local magnetization reads 



(S z ) =5-21n 



J 2 V5 
Ji 



regular terms. (12) 



The correction to the classical value 5 is a remnant of 
the logarithmic divergence in one dimension that has 
been regularized by an infrared cutoff w m i n = 0(J\VS) 
brought about by the interchain coupling Hl()(l . The ar- 
gument of the logarithm is the ratio of the maximum 
magnon frequency w max = 0( J2S) to the infrared cutoff. 
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Taking this formula at face value we estimate that the 
Neel order is present if the interchain exchange exceeds 
a critical value 

•he = 0(J 2 VSe- s/2 ). (13) 

Note that the ordering is a truly collective effect since 
independent spin chains possess no LRO. This feature 
makes the order-by-disorder problem rather different 
from its higher-dimensional analogues where each indi- 
vidual unit (say, spin plane in a canonical example of 
two interpenetrating square lattices) is ordered even in 
the absence of frustrating inter- unit interactions. Finally, 
Eq. H13|) indicates that at a large S there is a narrow re- 
gion of ratios J1/J2 without Neel order. This result is in 
agreement with a weak-coupling renormalization group 
analysis of the 5 = 1/2 problem by one of us^ who ar- 
gued in favor of a gapless sliding Luttinger liquid ground 
state in a wide interval < J1/J2 < 0.8. 

III. CROSSED CHAINS: Ji < J 2 

The four vacua [Fig. I3b)-(e)] found in the limit of 
weakly coupled chains, J\ <C J2, remain the ground 
states of the system for all finite ratios J1/J2 < 1. To 
confirm this, we have computed numerically spin- wave 
spectra of all collinear classical ground states for a lat- 
tice 16 x 16 with periodic boundary conditions. Modulo 
the global 0(3) spin symmetry, there are 2 16 = 65536 
spin configurations to consider. The energy of zero-point 
motion J5J is indeed minimized by the four Neel states 
(|S1-||5J|. In fact, not only they remain degenerate (with nu- 
merical accuracy): their spin- wave spectra are identical. 
That surely means that there is a hidden symmetry at 
work. 

The observed fourfold degeneracy is caused by a spe- 
cial gauge- like symmetry discovered by Henley^ It exists 
whenever a lattice can be split into corner-sharing units 
(tetrahedra on the pyrochlore lattice, second-neighbor 
pairs in the present case) with total zero spin in any 
ground state. Nonzero eigenfrequencies of such a sys- 
tem can be obtained by solving the equations of mo- 
tions for the transverse components of the total spins 
of these units (Recr a , Ima^), which have the following 
simple formi^Si^i 

* ® a ~ JupSapdp. (14) 

Here S a /3 is the ordered (longitudinal) component of the 
spin shared by units a and /3. It can now be seen that 
whenever two ground states are related by an Ising gauge 
transformation 

S'a0 = -^aSapAp , A a = ±1, (15) 

their nonzero modes are also related, 

<r' a =A a o- a , (16) 



and have identical frequency spectra. Therefore gauge- 
equivalent ground states have the same zero-point energy. 

In the current context, sites a of the dual lattice 
are "tetrahedra" (squares with crossings) of the original 
checkerboard lattice. For equal exchanges J2 = J\, a unit 
contains four spins whose total spin vanishes in a ground 
state. For strong diagonal chains (J2 > Ji), the total 
spin must vanish on both diagonal bonds separately, so 
that there are two units on every site of the dual lattice: 
the diagonal links f and r\. The resulting equations for a 
collinear Neel state read 

%V V = S T ,r+f]{JlLlr+fi + ^^r+77) 

+ S r , r -fj(Jill r -fj + J2^r- jj), (17) 

where r are coordinates of a tetrahedron. A spin labelled 
S rr / is shared by the tetrahedra located at r and r'. Fi- 
nally, /ji T and v T are transverse spin components of its 
units: 

S m _i )fl + S m+ i !n = (Re^ mn ,lm^ m „,0), 
S m ,„_i + S mn+ i = (Rei/ m „,lmtw,0) (18) 

in the notations of Fig. 

The transformation (|15H16|) does not actually reflect a 
local symmetry: applied to a single unit, a, it flips the 
spins S a {3 shared by a with other units /3. These other 
units acquire a nonzero total spin and violate the ground- 
state condition. Therefore the transformations must be 
made on a number of dual sites (an infinite one for an 
infinite lattice). For J 2 > J\, entire diagonal chains of 
spins must be flipped. It can be checked that the two 
ground states shown in Fig. [21 (b) and (d) are related 
through such a gauge transformation flipping spins on 
every other diagonal chain in both directions. 

In addition to breaking the spin 0(3) symmetry, the 4 
ground states also violate the spatial symmetry of the 
checkerboard lattice. This leads to interesting conse- 
quences. Although the spin symmetry must be restored 
at any finite temperature (the Mermin-Wagner theo- 
rem), the discrete lattice symmetries need not. In such 
a case, the low-temperature phase can have a long-range 
spin-Pcicrls (bond) order. An example of such behavior 
was discovered by Chandra et alM for the Heisenberg 
antifcrromagnet on the square lattice with large second- 
neighbor coupling. 

To see the pattern of bond order in the proposed spin- 
Peierls states one can look at the bond averages (S^ ■ Sj) 
in the ground states of Fig. [5] (b)-(e). We have shown 
the frustrated bonds (those with parallel spins) in color 
depending on the bond orientation: red for vertical and 
green for horizontal ones. The lower part shows the dual 
lattice of "tetrahedra" each painted in the correspond- 
ing color. Using this language, the ground states can 
be described in a simple manner: second neighbors on 
the dual lattice have the same color. In fact, if the dual 
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sublattice is divided into two sublattices, each sublattice 
exhibits "ferromagnetic" Ising order — in terms of these 
color variables — independently of the other sublattice. 
Hence a fourfold degeneracy mentioned above. 

An ordered state thus can be completely characterized 
by a composite order parameter Z2 x Zi x S2' two in- 
dependent Ising variables and a Neel vector. The Ising 
order parameter, 

h = ((S r>r+ | - S P) r _ f ) • (S r , r+i5 - S r , r ^)) , (19) 

has a counterpart in the frustrated antiferromagnet on 
the square lattice with a second-neighbor coupling. 34 
However, we will find a richer structure of ground states 
because our order parameter is, in fact, a component of 
a doublet f = (/1, /2), defined for every tetrahedron: 

((S x + S 2 )-(S 3 + S 4 ) - 2(Si • S 2 + S 3 ■ S 4 )) 

h ' 

((S 1 -S 2 )-(S 3 -S 4 )) ( . 

[Fig. [3(a)]. The other component, /1, comes into play 
when the vertical and horizontal bonds become (nearly) 
equivalent to the diagonal ones, a situation encountered 
on the three-dimensional pyrochlore latticed 

IV. PLANAR PYROCHLORE: ,h = J2 

This is a point with a very large classical 
degeneracyi^S^i Only the total spin of a "tetrahedron" 
must vanish in a ground state, but not necessarily the 
spins of second-neighbor pairs separately. To the next 
order, 0(1/ S), numerical comparison of zero-point fluc- 
tuation energies in collinear Neel states still reveals a 
large degeneracy — much larger than in the previously dis- 
cussed case J 2 > Ji- The ground states of Fig. |21 (b)- 
(e) become degenerate with the Neel state of the simple 
square lattice [Fig. |2Ja)] and many others numbering 2 L 
in total. Apparently this multitude of degenerate ground 
states kills the long-range Neel order. On the other hand, 
it will be seen that the bond order survives. 

To proceed, we present an explicit construction of 
all collinear ground states degenerate at the 1/S level, 
and identify a short-range interaction that selects these 
ground states. It turns out that the most economical 
description of these states is obtained in terms of the 
bond — rather than spin — variables. 

A. Gauge-equivalent collinear states 

For equal exchanges, Eq. I|14f) holds for quartets of 
spins on "tetrahedra" of the checkerboard lattice (squares 
with crossings). As far as collinear Neel states are con- 
cerned, there are now three distinct possibilities: parallel 
spins can be found on vertical, horizontal, and now also 
diagonal bonds, which we encode, respectively, as red, 



green, and blue states of a tetrahedron. In the Neel state 
of the simple square lattice [Fig.|2Ja)] all diagonal bonds 
have parallel spins, so that this state is uniformly blue. 
After casting Henley's gauge principle in bond language, 
we will readily reproduce all collinear Neel states degener- 
ate with the blue vacuum. The four ground states found 
in the J 2 > Ji case [Fig. |2 (b)-(e)] are among these. 

As before, a gauge transformation on a tetrahedron 
involves flipping all its spins. Parallel spins are found 
on the same bonds before and after the transforma- 
tion, therefore the color of that tetrahedron remains un- 
changed; it is the colors of its neighbors that are affected. 
Therefore, gauge transformations can be done separately 
and independently on the two sublattices of tetrahedra. 

Flipping the four spins S a p on tetrahedron a from sub- 
lattice A takes four adjacent tetrahedra (3 E B out of the 
ground state. To fix this problem, for each (3 we must 
perform at least one more gauge transformation on one 
of its neighbors 7 G A. Here are the rules for gauge trans- 
formations performed on sublattice A: 

(a) For every tetrahedron (3 £ B, the number 
of gauge transformations A a = —1 on adja- 
cent tetrahedra a € A can be 0, 2, or 4. 

(b) If this number is or 4, (3 remains blue. 

(c) For 2 gauge transformations, the two a 
cannot be on opposite sides of (3 (e.g., north- 
east and southwest). 

(d) If both a are north of (3 (or both are south 
of 0), (3 becomes red. If both a are east (west) 
of (3, it turns green. 

By using these rules, we can now construct an arbi- 
trary ground state starting with the blue one. As the 
ground state is a direct product of independent ground 
states on sublattices A and B, we will construct a ground 
state of sublattice B by making gauge transformations 
on sublattice A. Suppose, for definiteness, that there is a 
red tetrahedron on sublattice B and that the two gauge 
transformations were made north of it [Fig. Efa)]. If no 
other gauge transformations were made, this state would 
violate Rule (a): the two B tetrahedra shown in black 
are not in their ground states. Additional gauge trans- 
formations cannot be made around the original red site, 
for it will become blue [Fig. Hfb)], contrary to the initial 
assumption. The only remaining possibility is shown in 
Fig. life): additional gauge transformations are made on 
the same horizontal line. Continuing the process we find 
a line of red tetrahedra extending over the entire sub- 
lattice B in the horizontal direction [Fig. 0{d)]. Thus a 
generic ground state of one sublattice consists of horizon- 
tal red and blue stripes of arbitrary widths [Fig. E{e,f )] 
or of vertical green and blue stripes [Fig. E£g)] . 

The process is then repeated with the roles of the sub- 
lattices reversed: sublattice A is colored via gauge trans- 
formations on sublattice B (whose colors are unchanged) . 
A sample ground state is shown in Fig. 0^h) . 



Now it is easy to count the number of degenerate 
ground states on an L x L checkerboard lattice with peri- 
odic boundary conditions. On a single sublattice of tetra- 
hedra, there are 2 x 2 i / 2 ~ 1 red-and-blue ground states: 
the exponential reflects the number of ways to place hor- 
izontal domain walls separating red and blue stripes; the 
prefactor accounts for a duplicate set of states with red 
and blue domains exchanged. In addition, there is an 
equal number of green-and-blue states bringing the total 
to 2 L / 2+1 — 1 per sublattice. (The blue state has been 
counted twice.) The total degeneracy of the ground state 
(including both sublattices of tetrahedra) is therefore of 
order 2 L . 



B. Long-range bond order 

Is there a spontaneously broken symmetry? There's 
definitely no long-range vector order. It is easy to see 
that two spins located on the same diagonal can be par- 
allel and antiparallel with an equal probability: they are 
parallel if there is an even number of nonblue tetrahe- 
dra in between, antiparallel if the number is odd. By 
construction (Fig. 0} , the probabilities of these outcomes 
are equal and (S r • S r /) = for these two spins: there 
is not even a short-range order. The argument can be 
extended to (almost) any other direction. (Exceptions 
are the vertical and horizontal directions: there are long- 
range correlations of spins along the horizontal or vertical 
lines of bonds in the cyan and magenta states.) 

The only remnant of the spin order is the collinearity 
of spins: (S r • S r /) 2 = S 4 in any ground state for any pair 
of spins. Thus one can conclude that the ground state 
is a spin nematic whose order parameter is a director 

(syz 2 ). 

In addition, the system has long-range bond order. In- 
deed, if one tetrahedron on sublattice A is colored red, 
there are no green tetrahedra anywhere on the same 
sublattice — and vice versa. Thus the symmetry between 
red and green colors (the symmetry between vertical and 
horizontal bonds) is spontaneously broken. The average 
color of a given tetrahedron is magenta (red and blue 
stripes) or cyan (green and blue). Each sublattice finds 
itself in the cyan or magenta phase, Figs. ^i,j). In the 
cyan phase, the vertical ("red") bonds feature antipar- 
allel spins; in the magenta phase, antiparallel spins are 
found on horizontal ("green") bonds. The colors of the 
sublattices are independent at order l/S. 

This argument can be made more precise by turn- 
ing to the bond order parameter f 2 (|20|) . As shown in 
Fig. OJa), f 2 is positive, zero, and negative in the red, 
blue, and green states, respectively. Averaging over the 
entire ground state manifold gives (f 2 (r)) = for any 
tetrahedron. The reasoning put forward in the previous 
paragraph suggests that (/^(r)/^!"')) = S 4 /2 > even 
as |r — r'| — > oo (as long as the tetrahedra r and r' reside 
on the same sublattice). Thus a long-range bond order 




FIG. 4: Construction of ground states at the order l/S for 
J% = Jz. The secondary color magenta is an equal-part mix- 
ture of red and blue; similarly, cyan is an equally-weighted 
average of green and blue. The letter G marks the locations 
of Z2 gauge transformations. Filled triangles represent tetra- 
hedra not in a ground state. 



is present 
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The long-range order is thus similar to the case J\ < 
J 2 , with one exception: the Neel vector in the composite 
order parameter is replaced with a director that defines 
a collinearity axis. The order parameter now has the 
structure Z2 x Z2 x S 2 j Z^- two Ising order parameters 
in addition to a spin nematic. The spin-nematic order 
breaks a continuous rotational symmetry and therefore 
will be lost at any finite temperature. The remaining 
Ising orders Z2 x Z2 are expected to survive up to a 
finite temperature 0{JS). 



C. Effective spin interaction 

One may wonder what kind of a Hamiltonian gives 
rise to a strongly degenerate set of ground states de- 
scribed above. In fact, it can be derived following Hen- 
ley's method. 1 - The quantum correction of order l/S 
to the ground state energy (|2J) is obtained by rewriting 
Eq. 114(1 for the eigenmodes: 

h 2 u 2 <7 a = J 2 22 E S a pSp 1 <7 1 

0(a) 7 03) 

= AJ 2 S 2 a a + J 2 J2 Yl S a pS M a^{2\) 

/3(a) 7 (/9)A» 

The notation (3(a) indicates that tetrahedron (3 is a 
neighbor of tetrahedron a (they share spin S a p)- Fol- 
lowing Henley, we introduce an adjacency matrix, 



^t»7 — 25*2 E SafiSp-y, 
18(0,7) 



(22) 



whose matrix elements are nonzero when tetrahedra a 
and 7 have a common neighbor (3. The sum over nonzero 
eigenfrequences can now be expressed in terms of the 
matrix T: 



E 



H\lu 



2 =JSTrVT^T = JSTrll--- T - 



T T 



This expansion converges rather slowly (eigenvalues of 
T extend all the way up to 1) and cannot be used for 
quantitative purposes. Nevertheless, it provides correct 
qualitative answers, as will be seen shortly. 

Because T aa = 0, the first spin-dependent correction 
comes from the second-order term, 



JS 



TrT 



-3 JSL 2 



— y 

apyo 



Sa/lSpjSjsSsa- (24) 



(In the last sum, all tetrahedra are distinct.) This spin- 
dependent term can be considered as a four-body inter- 
action of spins or as a two-body interaction of bond vari- 
ables. Indeed, the interaction —SapSp^S^sSsa is mini- 
mized when the number of down spins is even, which can 
be cast in bond language as 



, JSL 



FIG. 5: Ground states of Henley's four-spin interaction, 
Eq. (1241 . Frustrated bonds are shown in color. Parametriza- 
tion in terms of bond (rather than spin) variables decouples 
tetrahedra of one sublattice (a, 7) from those of the other 
(13,8). 



i.e., when the vertical bonds of tetrahedra (3 and 5 si- 
multaneously have parallel or antiparallel spins. All such 
bond configurations are shown in the top row of Fig. 03 
Alternatively, this 4-spin term can be considered as a 
bond interaction between tetrahedra a and 7, which is 
minimized when 



SfictS a p — SpjS-yS — iS" 



(26) 



i.e., when their horizontal bonds simultaneosuly have par- 
allel or antiparallel spins. These ground states are shown 
in the bottom part of Fig. |S| 

Note that the bond interaction generated by the term 
—TrT 2 is between tetrahedra of the same sublattice. 
Higher-order spin loops —TrT™ can also be represented 
as an interaction between n tetrahedra of the same sub- 
lattice. In terms of bond variables, there is no coupling 
between different sublattices at the order l/S. 

Before deriving the three-body interaction, it makes 
sense to check for the ground states of the two-body bond 
Hamiltonian that we have just obtained: the energy is 
— 1 for states shown in Fig. [S] and +1 for the remaining 
states. More concisely, 

For tetrahedra on the same sublattice, near- 
est neighbors in the horizontal direction need 
to be both red or neither red; in the vertical 
direction both green or neither green. 

It is easy to see that this rule gives precisely all the 
ground states found in the beginning of this Section by 
means of Henley's gauge argument — and no other states. 

It is quite remarkable that a crude two-body approxi- 
mation, 



TrVl - T « Tr(l - T/2 - T 2 /8), 



(27) 



SapSpj — SyfiSsa — iiS 1 



(25) 



correctly reproduces all the ground states. This fact 
seems to indicate that the neglected n-body interactions 
(which are by no means small) can be expressed in terms 
of the two-body potentials identified above. If this is the 
case, the many-body interactions shift all ground states 
by the same amount without breaking their degeneracy. 
(We have checked that the 3 and 4-body interactions are 
indeed reducible.) 



D. Effective bond interaction 

The four-spin interaction l|24|) can also be written in 
terms of the bond variables (|20|l . Adjacent tetrahedra a 
and 7 of the same sublattice (Fig. interact with energy 



F - JS ( l 



S 2 V3 



£7 * 

S 2 V3 



The unit vector n a7 points in the red direction [Fig.lsja)] 
for neighbors along x and in the green direction for neigh- 
bors along y. The one-body piece can be viewed as a 
"magnetic" field for the two-component "spin" f a that 
points in the cyan or magenta direction. Once we sum 
over all neighbors of a given tetrahedron, the resulting 
"magnetic field" has the blue color and the magnitude 
J/4S 3 \/3. The two-body potential couples ferromagneti- 
cally the red or green components of neighboring "spins" 
depending on the direction. 

It is worth noting that a very similar interaction 
has been obtained for the quantum case S — 1/2 by 
Tsunetsugui^Siii We will return to this point later in 
Sec. EH 



V. LARGE S: A BRIEF SUMMARY 

Our study of the ground state of the Heisenberg 
checkerboard antiferromagnet to order 1/S establishes 
the existence of a long-range bond order in this system. 
The bond order breaks a mirror symmetry of the lat- 
tice exchanging the vertical and horizontal directions (the 
green and red colors in our notation). The discrete (Zi) 
nature of the broken symmetry assures survival of the 
long-range bond order to finite temperatures. It is in- 
structive to trace the evolution of the ground state as 
the ratio of the first and second-neighbor exchange cou- 
plings varies. 

For J2/J1 < 1, the ground state is unique (up to a 
global rotation of all spins): it is the Neel state of the 
simple square lattice (J2 = 0). Although the global 0(3) 
spin symmetry is broken at T = 0, it is restored at any 
finite temperature. The Z2 symmetry is manifest. 

For J2/J1 > 1, we find a fourfold degeneracy as each 
sublattice independently chooses one of the two ground 
states with collinear spins. In the red state, parallel spins 
are found on vertical bonds (and antiparallel spins on all 
other bonds). In the green state, parallel spins are found 
on horizontal bonds (Fig. |2J). 

The point J2/J1 — 1 is special: there are 0(2 L ) de- 
generate ground states on an L x L lattice to order 1/S. 
There is no Neel order even at zero temperature. Still, 
long-range bond order is present: if a single tetrahedron 
is put in the red state, there are no green tetrahedra on 
the same sublattice. By averaging over all red-and-blue 
states (Fig. 0} one obtains a magenta phase in which 
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tal bonds. (For the rest of first and second neighbors, 



(b) 



FIG. 6: (a) Ground states computed to order 1/S as a func- 
tion of J2/J1. The two sublattices of tetrahedra are colored 
indepe ndently of e ach other, (b) Ground states for S — 1/2 
(Refs. 19 25 26 27). The secondary colors (cyan, magenta, 
and yellow) encode the location of satisfied bonds (those with 
antiparallel spins). They are the opposites of the primary- 
color states (respectively, red, green, and blue). 



(Si ■ Sj) — 0.) The cyan phase, where (Si • Sj) — —S 2 on 
vertical bonds and on the rest, is obtained by averaging 
over the green-and-blue states. 

The independent ordering of the two sublattices of 
tetrahedra is probably accidental: the quantum correc- 
tion of order 1/S couples bond variables of the same 
sublattice. Higher-order corrections may well produce 
intersublattice couplings locking the sublattice order pa- 
rameters. 

As we have shown for J2/J1 — 1, parametrization of 
the four-spin interaction H24JI in color terms gives a Potts- 
like model with q = 3 states and direction-dependent 
interactions (Fig. [SJ) . Small deviations of the difference 
J2/J1 from 1 amount to adding a magnetic field h pro- 
moting or suppressing the blue state. For J2/J1 < 1 
(blue field), it is advantageous to place parallel spins on 
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the weaker diagonal bonds making the entire lattice blue. 
For J2/J1 > 1, the field points in the opposite — yellow — 
direction suppressing the blue state and forcing the sys- 
tem to choose between the red and green states. Lastly, 
J2/J1 = 00 is the regime of decoupled diagonal chains. 
The phase diagram is shown schematically in Fig.EIa). 



VI. QUANTUM LIMIT: S = 1/2 

It is interesting to compare our large-5 answers to the 
numerical results for 5 = 1/2 obtained recently i 25 ' 26 ' 27 
For Ji = J 2 they have found a bond-ordered state in 
which the probability of finding a spin singlet is enhanced 
on half of the squares without crossings, the "quadrumer" 
phase in Fig. Efb). A similar plaquette state has ap- 
peared in the analyses of the quantum dimer model 
for the planar pyrochloreS and the Shastry-Sutherland 
latticed In our terminology, this plaquette state cor- 
responds to the magenta-cyan vacuum [Fig. EJa)]: the 
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bonds for one sublattice of tetrahedra and on horizontal 
bonds for the other. The 5=1/2 ground state quantum 
wavefunction, however, cannot be expressed as a simple 
product of single tetrahedron configurations. 

Although the large-5 and 5 = 1/2 answers are the 
same at equal couplings, there are also important differ- 
ences that become apparent when we compare the results 
for Ji ^ J 2 . In the large-5 phase diagram [Fig.HJIa)], the 
magenta-cyan state is a single point sandwiched between 
the blue and red-and-green phases. For 5 — 1/2, it is 
found in a finite range of ratios J2/J1 around 1. The 
phase diagram for 5 = 1/2 inferred from the works of 
Lhuillier et al. is shown in Fig.^b). 

The 5 = 1/2 case does not lend itself to a straightfor- 
ward analytical treatment. Nevertheless, we have found a 
useful phenomenological approach that sheds some light 
onto its phase diagram. The bond variables f , which we 
have introduced previously for classical spins l|20|) , can be 
defined in the same way for any spin value 5. The ground 
states of an isolated tetrahedron with a Heisenberg in- 
teraction between its four spins are 25 + 1 degenerate 
singlets. The operators f = (fi,f 2 ), 
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FIG. 7: Transitions between the Neel and "quadrumer" 
states, (a) Spin order melts at J2 / Jr = j a , bond order appears 
at J2/J1 = jb > js- (b) Spin order becomes noncollinear at 
J2/J1 = jb, then disappears at J2/J1 = j s > jb- 



would be proportional to the scalar product f Q -f^. In the 
5=1/2 case, 



fa 'f/3 = ~^J ax Tp x + T az T0 z ). 



(31) 



fx = 



(Si + S 2 )-(S 3 + S 4 ) - 2(Si • S 2 + S 3 • S 4 
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(S 1 -S 2 )-(S 3 -S 4 ) 



(29) 



leave this manifold of states invariant. In this Hilbcrt 
space they act as Hermitian matrices (25+1) x (25+l)4i 
In particular, for 5=1/2 they are proportional to the 
Pauli matrices. One can choose a basis in which 



h = —r x 



f V3 

h = —T Z 



(30) 



If the interaction between the bond variables on differ- 
ent tetrahedra a and f3 were of the pure Potts form, it 



However, as we have seen in Sec. IIV CI the interaction 
has a more complicated form. It is direction-dependent 
and is asymmetric in the Potts flavors: red and green are 
different from the blue even when J\ = J 2- The crudest 
way to reflect this asymmetry is to write an interaction 
of the form K 1 f al f /31 + K 2 f a2 f/32 with K x ^ K 2 . In 
addition, there can be a "magnetic field" coupling —hf a i 
that selects the blue states (for J 2 <C Ji) or the red and 
green states (J 2 3> Ji) in the classical case. 

H = ^2 [Kxx{Ra0)TaxT0x+K zz (R al3 )T az Tp z }—h'^2 T ax- 

(a/3) a 

(32) 

Here we assume for simplicity that the potentials K xx 
and K zz depend on the distance between tetrahedra, but 



11 



not on the direction. Our findings in the previous sec- 
tion indicate that the second-neighbor coupling is ferro- 
magnetic (K < 0). Numerical results of the Fouet et 
alM> suggest that the nearest-neighbor interaction is an- 
tiferromagnetic (K > 0). The fact that the Ising order 
parameter is always related to the component f% points 
to the dominance of the K zz coupling over K xx . In this 
case we can neglect the K xx terms altogether without 
changing the critical properties in any significant way: 

H = ^ K{R a p)T az T fjz - h^r ax . (33) 

(a/3) a 

This is the Hamiltonian of an Ising antiferromagnet in 
a transverse magnetic field, whose properties are well 
known. At h = 0, it is a classical system with long- 
range order: (r z ) = +1 on one sublattice and —1 on the 
other. (The tetrahedra form a square lattice, so that the 
interactions are not frustrated.) The Neel phase is pre- 
served in a finite range of weak magnetic fields \h\ < K . 
A strong "magnetic field" h 3> K induces fast quantum 
fluctuations that kill the Neel order even at zero tem- 
perature, so that (r 2 ) = 0. The zero-temperature phase 
diagram is shown in Fig. HJ1 (b). The system 1MMI) has 
an "antifcrromagnctic" phase for \h\ < h c , which corre- 
sponds to the "quadrumer" state of Fouet et al. 25 Quan- 
tum critical points at h — ±h c separate this phase from 
the quantum disordered regions with no bond order — the 
yellow and blue states in which the diagonal bonds are 
strong ((r x ) > 0) and weak (t x ) < 0), respectively. 

Let us pause for a moment and take a critical look at 
this phenomenology From the viewpoint of bond vari- 
ables, at some value J2/J1 = j c ~ 1, which corresponds 
to h = 0, conditions are ideal for a spin-Peierls state: 
quantum fluctuations induced by the "transverse field" 
h oc J 2 / J\ — j c are absent, so that bond order is robust. 
Going away from that point increases quantum fluctu- 
ations of f2 oc t z and reduces the bond order parame- 
ter until the bond order completely melts at a quantum 
critical point. If the critical behavior is adequately de- 
scribed by the Ising model, the energy gap vanishes at 
the critical point only. Both phases — bond-ordered and 
bond-disordered — are gapped. 

This may or may not be the case. Consider in more 
detail the transition from the Neel phase of the simple 
square lattice (J2/J1 = 0) to the spin-Peierls phase at 
J2/J1 ~ 1. There are three distinct possibilities: (1) As 
J2/J1 increases, first the spin order melts at J2/J1 = j s , 
then a bond order appears at J 2 / Ji = jb > j s , see Fig. \7\ 
(a). The intermediate phase has neither spin, nor bond 
order. The critical behavior near J2/J1 = jb is ade- 
quately described by our bond phenomenology. (2) Mag- 
netic order persists into the spin-Peierls phase. The spin 
and bond orders coexist in the range jb < J2/J1 < js- 
In this case, the emergence of bond order at J2/J1 = jb is 
more adequately described as a spin canting transition, 
see Fig. [7| (b) . The system is gapless on both sides of 
this quantum phase transition. (3) The spin order dis- 
appears simultaneously with the onset of the bond order 



at J2/J1 = j s = jb, as shown in Fig. (b). One phase is 
gapped, the other is gapless. 

As is the usual limitation of exact diagonalisations, the 
numerical data of Sindzingre et alm& were obtainable only 
for moderately small lattices (up to 6 x 6). On the basis 
of their data we are unable to tell whether or not the spin 
and bond orders overlap. We feel, however, that all three 
possibilities can be realised as a matter of principle, so 
the Ising phenomenology may be useful. 



VII. CONCLUSION 

In this paper, we have studied the ground state 
of a Heisenberg antiferromagnet with large spins S 
on the checkerboard lattice, also known as the pla- 
nar pyrochlorc. To zeroth order in 1/S — the classical 
approximation — the magnet has an extremely large, con- 
tinuous degeneracy of the ground state. In the next or- 
der in 1/S, this accidental degeneracy is partially lifted 
by quantum fluctuations. The main achievement of this 
work is a complete characterization of the ground-state 
properties of this magnet to order 1/S. 

The ground states with the lowest energy of zero-point 
motion are found among the classical vacua with collincar 
spins. By using a special dynamical symmetry discovered 
by Henley^ we have explicitly constructed all of these 
ground states and shown that their number is of order 2 L 
in a lattice L x L. We have shown that there is no long- 
range Neel order: in the ensemble of these ground states 
spin correlations (S^ • Sj) vanish beyond nearest neigh- 
bors. However, there is a long-range bond order: "tetra- 
hedra" of the checkerboard lattice spontaneously pick up 
one of the two states, in which the nearest-neighbor cor- 
relations (Sj-Sj) are uniformly — S 2 for vertical bonds 
and for horizontal ones — or vice versa. The bond order 
breaks the rotational symmetry of the lattice. 

More precisely, the bond order parameter is Z2 x Z2 : 
the two sublattices of tetrahedra order independently of 
each other, which is reminiscent of the Ashkin-Teller 
model. The discrete character of the broken symme- 
try indicates that the bond order will likely survive at 
low temperatures. The critical properties of the thermal 
transition between the bond-ordered and paramagnetic 
states remain an open question. 

In addition to the bond order, there is a nematic long- 
range order: every spin points along a common direction 
n or its opposite — n. 

Our large- S analysis is in reasonable agreement with 
numerical results for the S =1/2 planar pyrochlore anti- 
ferromagnet obtained by Fouet et alw&> and Berg et al£L 
They find a bond-ordered ground state of almost the 
same kind as we do. The ordering of the two sublattices 
of tetrahedra is no longer independent: opposite pat- 
terns of the bond order are chosen. The ground state at 
5 = 1/2 appears to be a doubly degenerate spin singlet. 
Numerical data suggest the presence of a spin gap, which 
seems to rule out the presence of spin order of the Neel 
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or nematic type. This brings up the question of stability 
of the nematic order that we have found at large S. A re- 
cent calculation of Canalsr^ based on the Dyson-Maleev 
approximation, indicates that collincar ground states are 
locally stable at large S but could become unstable be- 
low some critical value S c - Although the Dyson-Maleev 
scheme is not a controlled approximation for small S 1 , 
Canals' scenario is consist ent with the results reported 
here and in Refs. I25l26l27l 

A lack of symmetry between the first and second- 
neighbor bonds compells one to study a more general 
system with unequal first and second-neighbor exchanges 
Ji ^ J%- For large S, the deviation of j = J2/J1 from the 
critical value j c = 1 plays the role of a "magnetic field" in 
the 3-state Potts model. The 3 flavors correspond to the 
3 collinear ground states of spins on a tetrahedron (mod- 
ulo a global rotation of the spins). They can be labelled 
by the location of frustrated bonds: diagonal, vertical, or 
horizontal. The "magnetic field" h cx J2/J1 — 1 prefers 
the "diagonal" state when J2 < J\, in which case the 
ground state is unique and the lattice symmetries are 
intact. When J2 > Ji, the "vertical" and "horizontal" 
states are favored, which leads to a spontaneous break- 
ing of the rotational symmetry of the lattice. Still, even 
at the critical point J2/J1 = 1 the symmetry between 
the three flavors is not restored: the diagonal bonds re- 
main different from the horizontal and vertical ones, so 
that the symmetry is still Z2, rather than S3. Accord- 
ingly, the order parameter is of the Ising, rather than the 
3-state Potts, model. Put another way, only one compo- 
nent of the Potts order parameter is used: that which is 
orthogonal to the direction of the "magnetic field" h. 

At large S, Neel order is present on both sides of the 
critical point J2/J1 = 1, but not at the critical point it- 
self. For Ji < J2 (diagonal chains with a weak frustrated 
coupling) the Neel order is induced by quantum fluctua- 
tions of spins. There are signs that, for a sufficiently weak 
interchain coupling, J\ < J\ c — 0{ J2 VSe~ s / 2 ), the Neel 
order may be destroyed. The fate of the bond order is 
unknown, although chances are that it is less susceptible 
to long-wavelength spin fluctuations. 

While it may be unreasonable to expect quantita- 
tive information about small spin values from a 1/S ex- 
pansion, we were tempted to make some general state- 
ments about the observed behavior of the S = 1/2 
system^iSSiSl It appears that there is some family re- 
semblance and that the zero-temperature phase diagrams 
for large S [Fig. EJa)] and for S = 1/2 [Fig. Efl>)] can 
be understood in similar terms, as far as bond order is 
concerned. We find it plausible that the bond operators 
(rerpesented for ,5 = 1/2 by 2x2 Pauli matrices) behave 
as spins of an Ising antiferromagnet in a transverse mag- 
netic field. The Ising order parameter is an expectation 
value of t z ; the "magnetic field" h couples to t x . In zero 
transverse field, there is a bond order of the antifcrro 
type, as observed by Fouet et al. A nonzero transverse 
field h cx J2/ J\— jc induces quantum fluctuations of the 
"spins" t z . At some critical value of the "field" the bond 



order melts and the ordered phase has a finite extent, in 
agreement with the numerical work of Sindzingre et alm& 
It is rather intriguing to find a valence-bond solid with- 
out Neel order in the limit of large S. This result is 
probably not unique to the planar pyrochlore and we in- 
tend to pursue this avenue further by studying other lat- 
tices. In particular, it would be interesting to solve a sim- 
ilar problem on the three-dimensional pyrochlore lattice 
and determine the resulting phase. In that case, the S3 
symmetry of the Potts states is intact and the outcome 
should be different. There are further open questions. 
Are there lattices that contain valence-bond liquids at 
large S 1 ? What happens at intermediate values of S 1 ? We 
look forward to finding out. 
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APPENDIX: ORDER FROM DISORDER IN THE 
CHAIN LIMIT 

Here we derive the results presented in Sec. [H] 

1. Notation 

The starting point J\ = corresponds to completely 
decoupled chains running along the diagonal directions 
£ and r). For technincal reasons, we will consider chains 
with an Ising anisotropy: 

#o = J2 £ S[Sf S? +1 + SfSf +1 ] + S?S? +1 (A.l) 

i 

Introduction of the anisotropy (0 < 6 < 1) helps stabi- 
lize the collinear ground state of a single chain: at the 
Heisenberg point ((5 = 1) Neel order along the chain is de- 
stroyed by quantum fluctuations. At the technical level, 
this is caused by a divergent 1/S correction to the stag- 
gered moment for 5 = l4» Then, strictly speaking, the 
initial assumption of the long-range Neel order on the 
chain breaks down. Therefore we compute the effects 
of the interchain coupling at 5 < 1 and then take the 
Heisenberg limit 5 — > 1. Unlike the staggered magnetiza- 
tion, the energy of zero-point motion J2} does not have 
an infrared divergence. This justifies approaching the 
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Heisenberg limit S — ► 1 from below4^ (Note, however, 
that the interchain coupling is always taken to be of the 
isotropic, Heisenberg kind.) 

At the Ising point, 6 = 0, the system has an extensive 
degeneracy: there are 2 L ground states. We parametrize 
them by introducing a single Ising variable ±1 for every 
chain. Then, for example, on the mth chain along £ — 
whose spins have coordinates r = (n + 1/2, m) — we have 
S* = (-l) n+1 s m S. We will use the Holstein-Primakoff 
transformation keeping the terms of orders S, S 1 ^ 2 and 
1. For a spin with S z > 0, 



S z = 

s+ = 
s~ = 



S — a) a 
V2S 



1 



\/25a t + O ( 



For spins with S z < we rotate the reference frame about 
the direction x through ir. Then for a spin located on the 
rn-th chain running along £ we obtain 

1 

(-l) n+1 s m 
(-l)" +1 Sl 

/ 



A) 



(A.2) 



V 



S — ala 



where r = (n + 1/2, m). 

The effective potential generated by quantum fluctua- 
tions is a function of L Neel vectors n m . Invariance under 
global spin rotations implies certain restrictions on the 
form of that function. For instance, a pairwise potential 
coupling chains with staggered magnetizations fii and ri2 
must be a function of the scalar product fii • ri2 . In addi- 
tion, symmetry of the lattice requires — for chains running 
in perpendicular directions — that the coupling be invari- 
ant under fii fii, so that it must depend on (fii -n 2 ) 2 . 
With this in mind, we will compute the lowest-order 
effect — O(Jf) — as a function of L(L — l)/2 angles be- 
tween staggered magnetizations 9 mn = arccos (n m • n n ). 
The pairwise nature of this potential allows us to tilt the 
spins uniformly on all r\ chains. It is convenient to choose 
fii as the z direction and let n.2 lie in the yz plane. 

For potentials coupling more than two chains, one must 
consider a general orientation of the staggered magneti- 
zations involved. However, because the leading term (0J 
already selects collinear states, in higher orders we will 
work with collinear configurations only. 

Thus for the nth rj chain — whose spins reside at r = 
(n, m + 1/2) — we perform an additional uniform rotation 
in the yz plane: 



SZ 



1 

(-l) m+1 i„cos6» (-1) TO+I i n sin0 | 
-(-l) ro+1 i„sin0 (-l) m +H n cos6 



( ^(a F + 4)\ 

V s 



a'a r 



(A.3) 



Naturally, the intra-chain Hamiltonian is not affected by 
these unitary rotations: 

#o = J2S [ 2a Ur + S(a r a r+i + ojo^ )] (A.4) 

r 

for the £ chain. After a Fourier transform, 

Ho = J2S ^2 [ 2a p a p + <5 cosp(a p a_ p + a p a}_ p ) , (A. 5) 
p 

where p is the lattice momentum along the chain direc- 
tion. Thermal Green's functions are given by 



(a t (u),p)a(-u,-p)) = U ° 2 ; lUJ 2 
(a(w,p)a(-uj,-p)) - 



(A.6) 
(A.7) 



uj^ + e p 

where w is a bosonic Matsubara frequency and 
cj = 2J 2 S, 

7 P = Scosp, (A. 8) 

The perturbation term in the Hamiltonian (for which 
we introduce no anisotropy) 

n — l/2,m + S„+i/2 .7)1 n,m — 1/2 1 ^n,m+l/2 ) 

(A.9) 

couples linear combinations of spins S„_ 1 / 2 , m + ^n+i/2,m 
and S„ TO _i/ 2 + Sn,m+i/2- Note that both linear combi- 
nations live on the same tetrahedron centered at (n, m). 
It is convenient to introduce variables £ nm and r\ nm rep- 
resenting transverse spin fluctuations on the respective 
diagonal links of the tetrahedron. For a link along £, 



= 2VS[ s m (-l) n+1 Q 

5m (-l)" +1 ^ 



(A.10) 



By direct comparison to (|A.2|I we obtain 

tx 

S// in 



y 

it m 



f \ f , \ 

^/f K- 1 ,m + a n-i,rn + a n+ ±,m + a «+I,m): 

— =(a' 1 — a„_ 1 
2V2 ™"5' m 2 



= 0. 



(A.ll) 



The longitudinal component £ z is of order S^ 1 ^ 2 and has 
therefore been dropped at the current level of approxi- 
mation. For r\ chains we similarly define r) a fields via 

1 , 1 1 > 





1 


Inm 




Inm 


i 






= 0. 



(A.12) 
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The sum of the spins along an r/ link is then expressed, 
according to l|A.3|) . as 

/ Vnm \ 

S„, m _i + S„, m+ 1 = 2VS (-l) m+1 t n cos 9 rjy m (k ,13) 

The interchain coupling (|A.9|I now reads 

V = 4J 1 Sj2(& m vZ m +e nm V y n m ), (A.14) 

where we have introduced a shorthand notation 
f y = f-lVN P v f\ v — (-l) m t cos# n y 

Snm V ± J °mS,nm> 'Inm — V '■J L n > - >Jn w 'Inm- 

(A. 15) 

To complete the preparation stage we work out Green's 
functions of the £ and r/ variables for uncoupled chains: 

<0)& m ^')> = ^G ab (r-r',n-n'IA.16) 
KmMm'iT')) = 5 m ,G ab (r-r',m-m'^.17) 

where indices a and b take on values x and y. The space- 
time Green's functions G ab (r, n) are easily obtained in 
terms of their Fourier transforms, 



which are given by the matrix 
G(u>,p) = 



-iwr+ipn 



(A.18) 



ujo(1— 7p) cos (p/2) — iaj sin (p/2) cos (p/2) 



iu; sin (p/2) cos (p/2) u;o(l+7p) sin 2 (p/2) 



(A.19) 

Partial Fourier transforms G ab (ui,n) are real and satisfy 
the following identities: 

G xx (cj,n) = (-l) n G yy (Lu,n) = G xx (u,-n), 

G xy (u;,n) = G yx (u;,n) = ~G xy (u;,-n), 

G xy (u;, n) = (-l) n+1 G xy (w, n). (A.20) 

The last line suggests that the off-diagonal components 
vanish for even distances n. At the Heisenberg point, the 
diagonal components vanish for odd distances. 



angle 8 comes through the y components — see Eq. (IA. 151 . 
Furthermore, terms containing two x and two y compo- 
nents vanish by symmetry: they have one factor of ,s 
only and staggered magnetizations s m are odd under re- 
flections in any line of spins along rj. Therefore the only 
possibility is four y components: 

lev f.v eV f.V \ _ (_-\\m+m'+n+n' 2 a 
\<,nm'lnm<,n'm''ln'm' / — I 1 J LUb ° 

x S m S m 't„i„/ (^nm^n'm') i^nrnVn'm') ' (A. 23) 

where we have used Gaussian statistics of the variables £ 
and rj. It is evident that lattice points (nm) and (n'mf) 
must belong to a £ chain and an rj chain simultaneously, 
so that they are the same point. This term therefore 
gives rise to a contact interaction of crossing chains. The 
correction to the energy of a ground state is found by 
taking the limit of zero temperature ((3 — > oo): 



E 



(2) 



(4Jiff) s 
2! 



2^ 



[G ra (w,O)] 2 £cos 2 



j — 



(A.24) 



where n m and n„ are the directions of staggered mag- 
netizations on chains running in the £ and rj directions, 
respectively; G(io,0) is the real-space Green's function 
G(oj, n) at distance n — 0. The numerical constant 1(5) 
is given here for the Ising (6 = 0) and Heisenberg (6 = 1) 
limits: 



m = \ 



2tt 



(A.25) 



Eq. I|A.24(I constitutes an order- by-disorder effect: 
collinear spin configurations (n m = ±n„) minimize the 
energy of quantum fluctuations. There are 2 L such 
collinear states in an L x L lattice with periodic bound- 
ary conditions. Their degeneracy is partially lifted at the 
order jf , as we discuss next. 



3. Order jf 



2. Order j\ 

The first nonvanishing correction to the free energy 
comes at the second order in J 2 and can be expressed as 
the second moment of the Eucledian action: 

**~li((jf *")')■ ,A - 21) 

The perturbation V given by Eq. I|A.14(1 . Its second mo- 
ment contains quartic averages 

l(P x n x A-P y fi y MP X n x 4-P v fi y \\ 
\\>nm Inm ~ Snm Inm J VSn'm' In' m' ~ Sn'm' In'm' / / ' 

(A.22) 

(Time variables are omitted for brevity.) Dependence on 
the staggered magnetizations s m and i„ and the tilting 



In the remainder of this section we will consider 
collinear vacua. Therefore the unitary transformation 
matrix in Eq. i|A.3|) becomes diagonal and l|A.13|) acquires 
the same form as (jA.lOf) . The ground-state energy will 
depend on the Ising variables s and t. The first non- 
trivial correction to the energy comes from the term of 
the fourth order in J%: 




(A.26) 



The dependence on staggered magnetizations s and t 
comes through the vertex ^ v fj v (|A.15|I . Symmetry argu- 
ments used in the preceding calculation show that s and 
^-dependent diagrams must contain exactly two vertices 
£ x rj x and two £ y fj v . We obtain 
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= - 2 x 6 x M|L / ^ E [(-!) fcs ™ W OGW(w, fc)G*-(w, -Z)G*>, -fc) 

+(-l)<Mn+« G xx (w, l)G xy (u>, k)G™{u, -l)G yx {u, —k) 
+ {-l) l+k s mSm+k t n t n+l G xy (uj, l)G yx (u, k)G xy (LJ, -l)G» x (w, -k)] . (A.27) 



Feynman diagrams contributing to Eg. (|A.27|I consist of 
a rectangular path formed by two £ chains and two r\ 
chains (Fig. |SJ. The factors 2x6 have a combinatorial 
origin. At this order, the interchain coupling J\ generates 
a two-chain interaction (between parallel chains) and a 
four-chain one (involving four crossing chains): 

E(4) = 2 J£IJ2J:Y: s m s m+k (A.28) 

2 ;>o fc>0 m n 
+A[ t n t n +i — Bki s 

Dimensionless couplings A\ > and Bui > arc 

A k = 2 7 (J 2 Sf f |%*>,fc)] 2 Z)] 2 , 

B kl = 2 8 (J 2 S) 3 J ^-[G x v(u,k)] 2 [G xy (u,l)] 2 .(A.29) 

They fall off quickly with the distances k and I. We have 
used the properties of the Green's functions l)A.20|) in 
deriving them. Below we discuss the Ising (6 — 0) and 
Heisenberg (S — ► 1) limits. 

a. Ising limit 

In the Ising limit, the magnons £ and r\ have an infi- 
nite mass and are unable to move far along the chain. 
Therefore G ab {u; 7 1) = for |/| > 1. As a result, the effec- 
tive potential depends on the nearest-neighbor products 
^m^m+i and t m £ m _|_i. 

-^Ising = 128 "jS" ^ ^ ( s mSm+l + tntn+l 
2 m n 

t„t„+i). (A.30) 

It is minimized by the ground states ©, in which stag- 
gered magnetizations have opposite signs on neighboring 
parallel chains. Such configurations minimize every term 
in Eq. (1A~30)| . 

b. Heisenberg limit 

At the Heisenberg point, diagonal (off-diagonal) com- 
ponents of the propagator (|A.19|) contain even (odd) har- 
monics of the translation operator e lk only. In real-space 




FIG. 8: Computation of the quantum correction to the free 
energy at O(Jf). Filled dots reprsent vertices £mnVmm open 
ones £mnVmn', the latter contribute a factor (— l) m+n s m in to 
the diagram, (a) The first term in Eq. 1A.27II . The diagram 
has a prefactor (— l) l t„t n +i. (b) The third term in Eq. 1A.27L 
Prefactor (~l) k+l s m s m+k t n t n+ i. 



terms, a magnon can travel an even distance by keep- 
ing its polarization, an odd distance by flipping it. Thus 
pairwise chain interactions are induced for parallel chains 
with an even separation; quartic ones involve pairs of par- 
allel chains with an odd separation: 

E(4) = 2J^^^^^ kSmSm+2k (A3i) 

2 k =l 1 = 1 m n 
+ A2lt m t m+ 2l — B-2 k ~\3.l-\ S m S m +2k-l tntn+2l-l)- 

The dimensionless coefficients A and B are given below. 
As previously noted, this potential is invariant under a re- 
versal of spins on every other diagonal chain JBJ. We will 
discuss the origin of this dynamical symmetry in Sec. IIIII 
and demonstrate that it holds to an arbitrary order in 
Ji. 

The dominant term in l|A.31|l — proportional to B\^\ — 
describes a coupling of four chains intersecting around 
an open square. It is minimized by spin configurations 
in which diagonal chains running in the same direction 
have either constant staggered magnetizations JSJ) , or al- 
ternating ones 10- These states — shown in Fig. [21 (b) 
through (e) — are related by a staggering transforma- 
tion JIjJi. Therefore they remain degenerate even upon 
the inclusion of all two- and four-chain interactions in 
Eq. IpOTJl . 

The second-largest term — a two-chain potential pro- 
portional to A2 — is minimized by Neel states of a dif- 
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ferent kind: the ones with opposite staggered magnetiza- 
tions for second-neighbor chains. Thus we obtain another 
viable candidate for a ground state: 



S4n 



Sin+1 



-SAn+2 



-Sin+3 



(A.32) 



and similarly for t. 

Because the fourth-order energy correction © contains 
oscillating terms, it is best to evaluate the energies at 
order J 4 exactly. We find that the energy of the states 
(JSJ and JJjJ is indeed lower by CJfS/ Jf per spin, where 



C = 



y/2 



45 
128 



1.991 x 10" 



(A.33) 



4. Coupling coefficients Ai and Bk 



In what follows we specialize to the Heisenberg case, 
5 = 1. Coefficients (|A.29ll of the two- and four-chains in- 
teraction (JSJ) are given in terms of the following integrals 



G xx (u>, n) 



uj r dp sin 2 pe ipn 



2 J_ 7T 2t: uj 2 + o)q sin 2 p ' 



2i 2n uj 2 + loq sm p 
Introduce an auxiliary variable u: 

uj = luq sinhu. (A. 35) 

We find: 



S n o — tanh \u\ 
G*>,n) = { 2o7 


where 5 mn is the Kronecker delta; 



e ! ntt ' n is even 



n is odd, 



n is even. 



(A.36) 



G xy (u,n) = { tanh it _,„„, , (A.37) 

' e |Tlii| sgnn n is odd. v ' 



2cjn 



tanh \ u\ 
8wq cosh 2 u 



(A.38) 



Substituting these into Eq. (|A.29|I produces the following 
expressions for nonvanishing coupling constants in the 



effective potential <|A.28|) : 

- 1 dx {l-xfx^- 1 ' 2 



Ai= I - 

Bki -- 



2^ 

1 dx (1_ 
2^ 



(1 



xp 

„|fe| + |i|-3/2 



I is even, (A. 39) 



k and I are odd. 



io m (l + x) 3 

Numerical values for the first few terms are: 

A 2 = 1.41 x lfr 3 , Bi A = 6.86 x 1(T 3 , 

A± = 1.63 x 1(T 4 , Bi <3 = 3.47 x 1(T 4 , (A.40) 

A 6 = 3.84 x 1CT 5 , B ltS = 5 3 ,3 = 5.29 x 1CT 5 . 
Although these coupling constants fall off rather 
quickly — as \l\~ 4 and (\k\ + \l\)~ 5 at large distances — 
partial cancellations in Eq. I|A.28|) call for a careful com- 
parison of energies of the candidate ground states. In a 
state with equal staggered magnetizations on all chains 
(JSJ, all products s m s m+ k = t m t m+ i = 1. The fourth- 
order correction to the energy is 

oo oo \ 




k=l 1=1 



(A.41) 

per spin. The state with alternating staggered magneti- 
zations © is degenerate with it by virtue of the stagger- 
ing symmetry ijfjj. In the (+H ) state <|A.32|) . s m s m+ k 

oscillates as a function of m (and thus averages out to 
zero) for an odd k and equals (— l) fe / 2 for an even k. 
Hence the expression for the correction to the ground- 
state energy at this order: 



E 



(4) 



N„ 



4*x2£>(-l)' 
2 i=i 



(A.42) 



The energy difference is then 



E 



- E 



(4) 



j4 Q 



(A.43) 



v spins 

where the constant is 

- 1 dx (1 - xf (1 - 4x + x 2 ) x 1 / 2 



C 



o 2tt 
V2 45 
~T ~ 128 



(1 +xf (1+x 2 ) 
s 1.991 x 1CT 3 . 



(A.44) 
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